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Abstract Astrophysical fluids, including interstellar and in- 
terplanetary medium, are magnetized and turbulent. Their 
appearance, evolution, and overall properties are determined 
by the magnetic turbulence that stirs it. We argue that ex- 
amining magnetic turbulence at a fundamental level is vital 
to understanding many processes. A point that frequently 
escapes the attention of researchers is that magnetic turbu- 
lence cannot be confidently understood only using "brute 
force" numerical approaches. In this review we illustrate this 
point on a number of examples, including intermittent heat- 
ing of plasma by turbulence, interactions of turbulence with 
cosmic rays and effects of turbulence on the rate of mag- 
netic reconnection. We show that the properties of magnetic 
turbulence may vary considerably in various environments, 
e.g. imbalanced turbulence in solar wind differs from bal- 
anced turbulence and both of these differ from turbulence 
in partially ionized gas. Appealing for the necessity of more 
observational data on magnetic fields, we discuss a possibil- 
ity of studying interplanetary turbulence using alignment of 
Sodium atoms in the tail of comets. 

Keywords Turbulence, MHD, Interstellar Medium 

PACS 95.30.Qd, 52.30.Cv, 96.50.Tf 



1 Introduction 

It is well known that astrophysical fluids are magnetized 
and turbulent (see Armstrong & Woo 1981, Armstrong et 
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al. 1995, Verdini & Velli 2007 and ref. therein). Turbulence 
is known to affect most of the properties of fluids, e.g. ther- 
mal conductivity, propagation of waves and energetic par- 
ticles, magnetic field generation etc. (see Kazantsev 1968, 
Moffatt 1978, Dmitruk et al. 2001, ScMickeiser 2003, Vish- 
niac, Lazarian & Cho 2003, Cranmer & van Ballegooijen 
2005 and references therein). Local Bubble, solar wind and 
interstellar medium are not exceptions: the fluids there are 
turbulent. A substantial progress in understanding of the me- 
dia above has been achieved using numerical simulations. It 
is very encouraging that present codes can produce simula- 
tions that resemble observations. However, one may wonder 
to what extent numerical results reflect reality. The answer 
to this question is not as simple as it may look. 

It is easy to justify a cautious approach to the interpreta- 
tion of the numerical simulations. Indeed, a meaningful nu- 
merical representation of the turbulent fluid requires some 
basic non-dimensional combinations of the physical param- 
eters of the simulation to be similar to those of the real ISM. 
One such is the "Reynolds number". Re, the ratio of the eddy 
turnover time of a parcel of gas to the time required for vis- 
cous forces to slow it appreciably. A similar parameter, the 
"magnetic Reynolds number", Rm, is the ratio of the mag- 
netic field decay time to the eddy turnover time. The proper- 
ties of flows on all scales depend on Re and Rm. It is not real- 
istic to expect that we can in any foreseeable future simulate 
turbulent flows with Re > 10^ and Rm > lO'^. Note that 3D 
simulations for 512 cubes can have Re and Rm up to ^ 2000 
and are limited by their grid sizes. It should be kept in mind 
that while low-resolution observations show true large-scale 
features, low-resolution numerics may potentially produce 
an incorrect physical picture. 

How feasible is it, then, to strive to understand the com- 
plex microphysics of astrophysical MHD turbulence? Sub- 
stantial progress in this direction is possible by means of 
"scaling laws", or analytical relations between non-dimensional 
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combinations of physical quantities that allow a prediction 
of the motions over a wide range of Re and Rm. Even with 
its limited resolution, numerics is a great tool to test scaling 
laws. 

On a basis of several selected examples we show that 
a "brute force" numerical approach has limited applicabil- 
ity while dealing with interstellar medium. Local Bubble 
and solar wind. Instead we claim that to create an appro- 
priate "tool box" that can be incorporated into simulations 
one must understand better underlying fundamental physi- 
cal processes. Moreover, we can define the circumstances 
when one should and should not rely on the present-day 
numerical simulations only by better understanding these 
processes. For instance, if magnetic reconnection is slow 
in collisional astrophysical environments, numerical simu- 
lations cannot represent fluids there. We argue below that 
in turbulent fluids reconnection is fast and therefore numeri- 
cal simulations are not very different from the reality in this 
respect. We believe that the solid progress in the studies of 
complex processes in magnetized astrophysical fluids must 
be based on verification of theoretical constructions with ob- 
servations. For this purpose we discuss a new technique of 
studying interplanetary magnetic fields using a subtle phe- 
nomenon of atomic alignment. 

We realize that within a relatively short review, one can- 
not cover in depth the properties and numerous implica- 
tions of astrophysical turbulence. Therefore while covering 
a few relevant subjects within the area our expertise, we pro- 
vide, wherever appropriate, references to extensive reviews 
and monographs. Although our review is intended for the 
special volume of the proceedings on the heliosphere and 
the Local Bubble, a number of examples in the review are 
from the area of interstellar medium. In fact, we feel that 
the existence of the gap between the interstellar and the he- 
liospheric/interplanetary communities is not healthy, as the 
fundamental underlying physics that researchers deal with 
in these two cases is very similar Moreover, we believe 
that with the Voyager spacecrafts approaching the interstel- 
lar medium (see Stone et al 2008, http://voyager.jpl.nasa.gov/) 
, the cross-polination between the fields should be enhanced. 

In what follows we discuss the spectrum and anisotropy 
of magnetic turbulence in different regimes, including the 
regime of balanced turbulence, imbalanced turbulence and 
viscosity-damped turbulence (§2). We consider turbulence 
intermittency in §3 and the modification of the magnetic tur- 
bulence in the presence of cosmic rays in §4. §5 is devoted 
to the selected implications of turbulence in interstellar and 
interplanetary medium, including scattering of cosmic rays, 
turbulent reconnection and perpendicular diffusion of cos- 
mic rays and heat. We present an example of interplanetary 
turbulence studies using observations of the sodium tail of a 
comet in §6. Our main results are summarized in §7. 



2 Spectrum of turbulence and its anisotropy 

2.1 Anisotropic MHD turbulence 

Magnetized turbulence is a tough and complex problem with 
many excellent monographs and reviews devoted to differ- 
ent aspect of it (see Biskamp 2003 and references therein). 
A broad outlook on the astrophysical implications of the tur- 
bulence can be found in a review by Elmegreen & Scalo 
(2004), while the effects of turbulence on molecular clouds 
and star formation are reviewed in McKee & Ostriker (2007) 
and Balesteros-Paredes et al. (2006). However, the issues of 
turbulence spectrum and its anisotropics, we feel, are fre- 
quently given less attention than they deserve. 

While turbulence is an extremely complex chaotic non- 
linear phenomenon, it allows for a remarkably simple statis- 
tical description (see Biskamp 2003). If the injections and 
sinks of the energy are correctly identified, we can describe 
turbulence for arbitrary Re and Rm. The simplest descrip- 
tion of the complex spatial variations of any physical vari- 
able, X{r), is related to the amount of change of X between 
points separated by a chosen displacement 1, averaged over 
the entire volume of interest. Usually the result is given in 
terms of the Fourier transform of this average, with the dis- 
placement 1 being replaced by the wavenumber k parallel to 
1 and |k[ = l/l'l- Foi" example, for isotropic turbulence the 
kinetic energy spectrum, E{k)dk, characterizes how much 
energy resides at the interval k,k + dk. At some large scale 
L (i.e., small k), one expects to observe features reflecting 
energy injection. At small scales, energy dissipation should 
be seen. Between these two scales we expect to see a self- 
similar power-law scaling reflecting the process of non-linear 
energy transfer 

Thus, in spite of its complexity, the turbulent cascade is 
self-similar over its inertial range. The physical variables are 
proportional to simple powers of the eddy sizes over a large 
range of sizes, leading to scaling laws expressing the depen- 
dence of certain non-dimensional combinations of physical 
variables on the eddy size. Robust scaling relations can pre- 
dict turbulent properties on the whole range of scales, in- 
cluding those that no large-scale numerical simulation can 
hope to resolve. These scaling relations are extremely im- 
portant for obtaining an insight of processes on the small 
scales. 

The presence of a magnetic field makes MHD turbulence 
anisotropic (Montgomery & Turner 1981, Matthaeus et al. 
1983, Shebahn et al. 1983, Higdon 1984, Goldreich & Srid- 
har 1995, see Oughton, Dmitruk & Matthaeus 2003 for a 
review). The relative importance of hydrodynamic and mag- 
netic forces changes with scale, so the anisotropy of MHD 
turbulence does too. Many astrophysical results, e.g. the dy- 
namics of dust, scattering and acceleration of energetic par- 
ticles, thermal conduction, can be obtained if the turbulence 
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spectrum and its anisotropy are known. As we discuss be- 
low, additional important insight can be obtained if we know 
turbulence intermittency. 

Estimates of turbulence anisotropy obtained in relation 
to the observations of magnetic fluctuation of the outer he- 
liosphere and solar wind (see Zank & Matthaeus 1992 and 
references therein) provided, for an extended period of time, 
the only guidance for theoretical advances. This resulted 
in a picture of MHD turbulence consisting of 2D "reduced 
MHD" perturbations carrying approximately 80% of energy 
and "slab" Alfvenic waves carrying the remaining 20% of 
energy (see Matthaeus 2002 and references therein). In other 
words, in the suggested picture the MHD turbulence was 
presented by two anisotropic components, one having wave 
vectors mostly perpendicular to magnetic field (the 2D one), 
the other having them mostly parallel to magnetic field (the 
slab one). This model became a default one for many calcu- 
lations of the propagation of cosmic rays (see Bieber, Smith 
& Matthaeus 1988, Bieber, Matthaeus & Smith 1994). On 
the contrary, guided mostly by compressible MHD numeri- 
cal simulations, the interstellar community adopted a model 
of the MHD turbulence where the basic MHD modes, i.e. 
slow, fast and Alfven are well coupled together and effi- 
ciently dissipate energy in shocks (Stone et al. 1998, Mac 
Low 1999). Little cross-talk between the two communities 
did not stimulate the interdisciplinary debates on the nature 
of MHD turbulence, which was regretful, as the magneto- 
spheric community has the advantage of the in-situ space- 
craft measurements. 

In spite of the intrinsic limitations of the "brute force" 
approach, we feel that reliable results can be obtained nu- 
merically if the studies are focused on a particular property 
of turbulence in order to get a clear picture of the underly- 
ing physics occurring on small scales ("microphysics") that 
cannot be resolved in "global" interstellar simulation^. 

We feel that it were the "focused" numerical simula- 
tions that allowed the validation of the Goldreich & Shridhar 
(1995, henceforth GS95) model of MHD turbulence. Indeed, 
GS95 made predictions regarding relative motions parallel 
and perpendicular to B for Alfvenic turbulence. The model 
did not predict the generation of any "slab" modes and, in- 
stead of pure 2D Alfvenic modes, predicted that most of 
the Alfvenic energy is concentrated in the modes with a so- 
called "critical balance" between the parallel and perpendic- 
ular motions. The latter can be understood within intuitive 
picture where eddies mixing magnetic field perpendicular to 



By contrast, numerical simulations that deal with many physical 
conditions simultaneously cannot distinguish between the effects of 
different processes. Moreover, they inevitably have a more restricted 
interval of scales on which energy is injected by numerics, initial con- 
ditions, or boundary conditions. Their results are, therefore, difficult to 
interpret in physical terms. 



its locals direction induce Alfvenic waves with the period 
equal to the period of the eddy rotation. This results in the 
scale-dependent anisotropy of velocity and magnetic pertur- 
bations, with the anisotropy being larger for smaller eddies. 

The relations predicted in GS95 were confirmed numer- 
ically for incompressibl^ (Cho & Vishniac 2000, Maron & 
Goldreich 2001, Cho, Lazarian & Vishniac 2002a) and com- 
pressible MHD turbulenc^ (Cho & Lazarian 2002, 2003, 
henceforth CL02 and CL03). (see Cho, Lazarian & Vish- 
niac 2003 for a review). They are in good agreement with 
observed and inferred astrophysical spectra. A remarkable 
fact revealed in Cho, Lazarian & Vishniac (2002a) is that 
fluid motions perpendicular to B are identical to hydrody- 
namic motions. This provides an essential physical insight 
and explains why in some respects MHD turbulence and hy- 
drodynamic turbulence are similar, while in other respects 
they are different. 

GS95 provided theoretical arguments in favor of low 
coupling between fast and Alfven modes, and low impact of 
slow modes to Alfven modes (see also Lithwick & Goldre- 
ich 2001). This challenged the paradigm accepted by the in- 
terstellar community. While the decomposition of MHD per- 
turbations into fundamental MHD waves was widely used in 
the Uterature (see Dobrowolny et al. 1980) it was usually as- 
sumed that the Alfvenic waves exist and interact with other 
waves for many periods (see a discussion in Zweibel, Hi- 
etsch & Fan 2003). This is not the case of the GS95 model 
of turbulence, where the Alfven modes non-linearly decay 
within one wave period. This reduces the time of interaction 
and therefore the coupling. Interestingly enough, in GS95 
model, the Alfvenic modes can affect slow modes, but the 
opposite is not true. These results were successfully tested 
in CL02 and CL03. 

Insight into fundamental properties of interstellar turbu- 
lence has already paid dividends. For instance, theories of 
cosmic ray propagation and dynamics of interstellar grains 
have been revised in view of better understanding of MHD 
turbulenc^ (see Cho & Lazarian 2005 and references therein) 
Similarly, the issue of turbulent support of molecular clouds 



- The notion of the direction being local is critical. Small eddies are 
affected by magnetic field in their vicinity, rather than a global field. 
No universal scalings are possible to obtain in the frame of the mean 
magnetic field. 

As in any developing field, ongoing controversies and competing 
ideas exist on how to improve the GS95 scalings (see Boldyrev 2005, 
2006, Beresnyak & Lazarian 2006, Gogoberidze 2007). 

^ Some studies of MHD compressible turbulence, e.g. Vestuto et 
al. 2003 did not perform a decomposition of MHD perturbations into 
Alfven, slow and fast modes as it is done in CL02 and CL03. They did 
not use local system of reference for which the GS95 scaling is for- 
mulated. Therefore a direct comparison of their results with the GS95 
predictions is difficult. 

' Incidentally, the concept of the scale-dependent anisotropy in re- 
spect to the local direction of magnetic field is veiy important. This 
local magnetic field is what actually affects energetic particles. 
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requires revisions. Indeed, a widely accepted view that the 
rapid decay of turbulence is caused by the coupling of com- 
pressible and incompressible motions was disproved with 
numerical simulations in CL02. In fact, it was shown there 
that the rate of turbulent decay was shown to depend on the 
ways that the turbulence is driven (Cho, Lazarian & Vishniac 
2003a). Moreover, contrary to common beliefs, MHD mo- 
tions do not die away at the scale at which ambipolar diffu- 
sion damps hydrodynamic motions. As we discuss in §2.1 a 
new regime of viscosity-damped MHD turbulence emerges 
(Cho, Lazarian & Vishniac 2002b). 

Some of the relevant results are illustrated in Figur^T] 
Contrary to the common expectation, the modes exhibited 
nice scaling laws that allow further analytical and numerical 
applications. For instance, numerical studies in CL02 and in 
CL03 reveal that the Goldreich-Shridhar scalings are valid 
for the Alfvenic part of the turbulence cascade even in the 
highly compressible regime (see also Beresnyak & Lazarian 
2006). 

The density spectrum of MHD turbulence was studied 
in Beresnyak, Lazarian & Cho (2006) and Kowal, Lazar- 
ian & Beresnyak (2007). The spectrum was shown to be ap- 
proximately Kolmogorov for low Mach numbers, and get 
flatter as the Mach number increases. The compression of 
gas by shocks was identified as the reason for the spectrum 
flattening. At the same time, it was shown that the loga- 
rithms of density have a Kolmogorov spectrum and GS95 
anisotropics. 

The statistical decomposition of MHD turbulence into 
Alfven, slow and fast modes suggested in CL02 was suc- 
cessfully tested for slow modes of magnetically dominated 
plasmas in CL03. The statistical procedure was further im- 
proved in Kowal & Lazarian (2006), where wavelets rather 
than the Fourier transforms were used. There the study of 
spectra and anisotropy of velocity, magnetic field and den- 
sity of compressible MHD turbulence was performed and 
the advantages of the wavelets in comparison with the Fourier 
technique were revealed for studying turbulence with weak 
mean field. For turbulence with B„,ean ^ SB the obtained re- 
sults are consistent with the CL02 and CL03 studies. Nev- 
ertheless, it would be wrong to say that we have a com- 
plete understanding of the scaling of MHD modes and their 
interactions. First of all, one should distinguish weak and 
strong Alfvenic turbulence. The weak turbulence is essen- 
tially 2D, with the turbulent cascade creating more structure 
perpendicular to magnetic field as the turbulence cascades 
(see Gaultier et al. 2000). Such a cascade emerges when the 
driving of turbulence at the outer scale is weak, i.e. the in- 
jection velocity is much less than the Alfven velocity. Al- 
though the weak turbulence picture corresponds to the early 
representation of MHD turbulence (see discussion in §2.1), 
one should keep in mind that the strength of Alfvenic in- 
teractions increases with the decrease of the scale along the 



cascade. Therefore the Alfvenic turbulence gets eventually 
strong, while both the inertial range and the astrophysical 
utility of the weak Alfvenic cascade are limited. The inter- 
action of weak Alfvenic turbulence with fast modes has de- 
pendences on the angle between B and the wavevector, as 
was shown by Chandran (2005). In addition, we discuss in 
§4 that the fact that no slab waves appear in MHD simula- 
tions does not mean that such modes do not appear in re- 
alistic astrophysical circumstances, e.g. when cosmic rays 
are present. Moreover, we show in §5 that in many circum- 
stances, the way cascade is initiated and dissipates affect 
both the properties and implications of magnetic turbulence. 

2.2 Imbalanced Turbulence 

MHD turbulence in the presence of sources and sinks gets 
imbalanced, in the sense that the flow of energy in one di- 
rection is larger than the flow of energy in the opposite di- 
rection. Solar wind presents a vivid example of imbalanced 
turbulence with most waves near the Sun moving in the di- 
rection away from the Sun. 

While theories of balanced MHD turbulence enjoyed much 
attention, the theory of imbalanced turbulenc^ attracted less 
work, unfortunately (see Biskamp 2003 and references therein). 
The analytical results were obtained for weak imbalanced 
turbulence (Galtier et al. 2002, Lithwick & Goldreich 2003) 
and they are applicable in a rather narrow range of imbal- 
ance ratios. Some earlier simulations of strong imbalanced 
turbulence were limited to rather idealized set ups (Maron 
&. Goldreich 2001, Cho, Lazarian & Vishniac 2002a), i.e. 
for the initial state the results of the simulations of strong 
balanced turbulence were used, but the amplitudes of waves 
moving in one direction were reduced. This did not allow 
making definitive conclusions about properties of imbalanced 
turbulence. 

We think that the best experimental data on the imbal- 
anced regime is currently available from observations of so- 
lar wind turbulence (e.g., Horbury 1999). This data, col- 
lected by spacecrafts, is consistent with Kolmogorov —5 /3 
spectrum, but does not provide sufficient insight into the 
anisotropy with respect to the local magnetic field. The im- 
balanced turbulence is not a rare exception, on the contrary, 
such processes as preferential decay of a weaker wave and 
the escape of turbulent energy from regions that generate 
perturbations make imbalanced turbulence widely spread in 
various astrophysical circumstances. 

Attempts to construct the model of stationary strong im- 
balanced turbulence were done in Lithwick, Goldreich & 
Sridhar (2007), Beresnyak & Lazarian (2008b), Chandran 
(2008). Below we discuss only the model in Beresnyak & 

^ Another name for imbalanced turbulence is a turbulence with non- 
zero cross-helicity. 
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Fig. 1 Statistics of MHD Turbulence: (a). Scaling of compressible motions for plasma with magnetic pressure ten times the gas pressure. This 
regime is important for molecular clouds (from CL02). The velocity (solid line) and B (dashed line) spectra are plotted against k (= l/(eddy size)). 
They show well defined statistical properties that allow further fruitful applications, (b). The simulations of decaying turbulence show that the rate 
of the decay of the total energy is a strong function of the imbalance of the energy contained in waves initially moving in opposite directions. In 
the lowest curve, the waves have the same amplitude. The energy decays more slowly with imbalance, affecting how long turbulence can support 
molecular clouds and whether transfer of energy between clouds is efficient (Cho, Lazarian & Vishniac 2002b). (c). The evolution of kinetic 
energy in 3-D compressible turbulence, initially started with V • v = as if the gas were incompressible. The dashed curve shows the evolution 
of the V • V = motions. We see that Alfvenic turbulence creates only a marginal amount of compressible motions, suggesting that Alfvenic 
modes should evolve independently of the compressible cascade (CL02). (d). Magnetic fluctuations persist beyond the turbulent damping scale at 
large k, while hydrodynamic fluctuations damp out in partially ionized gas (Lazarian, Vishniac & Cho 2004). This viscosity-dominated regime of 
turbulence may dominate small scale structure of partially ionized gas. (e-f). Isocontours of equal correlation for Alfven and fast modes (CL02). (e) 
The Alfvenic motions are much more correlated along B than perpendicular to it. (f) In contrast, fast magnetosonic fluctuations show essentially 
circular (isotropic) isocontours of correlation. 



Lazarian (2008b) as this is the only model that agrees with 
numerical simulations performed so far In view of a big 
picture, this model can be viewed as an extention of GS95 
model into the imbalanced regime. 

While the classic formulation of the GS95 critical bal- 
ance, based on causality, is unable to describe consistently 
the imbalanced case, Beresnyak & Lazarian (2008b) pro- 
posed a new way to introduce the balance between parallel 
and perpendicular modes. The new conditions for the criti- 
cal balanced between the parallel and perpendicular modes 
of oppositely moving waves was obtained appealing to the 
process termed "propagation cascading". In the case of the 
balanced turbulence the "old" and "new" critical balance 
condition results in the same GS95 scaling. However, in the 



case of imbalanced turbulence the new formulation actually 
predicts smaller anisotropy for the stronger wave, which di- 
rectly contradicts old causal critical balance, but is consis- 
tent with simulations. 

We assume that waves have different anisotropics, i.e. 
the dependence of longitudinal scale A to transverse scale X 
is different for each kind of wave. This situation is presented 
in Fig. 1, where some arbitrary longitudinal scale A ^ corre- 
sponds to the two different transverse scales, Ai for weak 
wave and X2 for strong wave w+. A+ is a longitudi- 
nal scale of w+ wave having transverse scale Ai. In case 
of strong turbulence, we expect that at least the is being 
strongly cascaded by w+. In this case the most effective mix- 
ing of on scale Ai will be obtained through w+ motions 
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Fig. 2 Left upper: w+ wavepacket, produced by cascading by wavepacket is aligned with respect to wavepacket, but disaligned with 
respect to the local mean field on scale Ai , by the angle 9. Left lower: the longitudinal scale L of the wavepackets, as a function of their transverse 
scale, /; Middle: the power spectrum of energies for both waves in an imbalanced forced incompressive 320"' numerical simulation. Right: the L(/) 
dependence in the same simulation, the length scales are in the units of the cube size. 



that are on the same scale. The longitudinal scale for will 
be provided by causal critical balance, since its cascading is 
fast. 

The cascading of w+ is somewhat more complicated. 
Since the amplitude of is not large enough to provide 
strong perturbations in w+, the w+ will be perturbed weakly, 
and the cascading timescale will be diminished according to 
the "strength" of the w^, just like it does in weak turbulence. 
Moreover, now the eddies will be cascading w+ eddies 
with similar longitudinal scales, which is the generic feature 
of weak cascading. 

The perturbations provided by will have a transverse 
scale of Aj. In other words, the energy of will be trans- 
ferred between A2 and Ai. The longitudinal scale for cas- 
caded w+ will be determined by the propagation critical 
balance in the following way. The wavepackets of are 
strongly aligned to the mean field on scale X\ and there- 
fore they are randomly oriented with respect to the mean 
field at a larger scale A2. The RMS angle of wavevector of 

eddies with respect to mean field on A2 will be around 
Q sa 5b'^ {X2) / V j^. This slant of wavepackets will deter- 
mine the increase of ^|| for newly cascaded w+ packets at Ai 
(see Figur^. The new interpretation of critical balance in 
the strongly imbalanced case is that the k\\ of the weak wave 
increases due to the finite lifetime of the wave packet, while 
in the strong wave it increases due to the field wandering 
of the strong wave itself on larger scales. This effect does 
not contradict the exact MHD solution of the wave propa- 
gating in one direction, because it requires the oppositely 
propagating wave as an intermediary. The de-alignment of 
the cascaded strong wave is possible because the weak wave, 
acting as a cascading agent, is strongly aligned with the field 
lines on scale which is different (smaller) than the scale of 
the strong wave it is acting upon. 

While the model in Beresnyak & Lazarian (2008b) is 
consistent with numerical simulations, higher resolution sim- 
ulations are definitely required. Moreover, the three mod- 



els of strong imbalanced turbulence, namely, Lithwick et 
al. (2007), Beresnyak & Lazarian (2008b) and Chandran 
(2008), provide predictions to be compared with the solar 
wind observations. Naturally, additional processes should 
be accounted for in the stellar wind studies, e.g. paramet- 
ric instabilities, reflections from preexisting density fluctu- 
ations (see Leroy 1980, Roberts et al. 1987, Hobery 1999, 
Del Zanna et al. 2001). This calls for further studies of im- 
balanced turbulence in compressible magnetized fluids. 

2.3 Viscous Turbulence 

When the magnetic turbulence takes place in viscous, but 
well conducting gas, its properties differ from those described 
above. The effect of the neutral gas can act in a number 
respects as the proxy of fluid viscosit}0. Schekochihin et 
al.(2004) and Goldreich & Sridhar (2006) argued that the 
plasma viscosity parallel to magnetic field can act in the 
same way as the normal viscosity of unmagnetized fluids (cp 
Braginskii 1965). If this is true, the fluids with high mag- 
netic Prandtl number Pr„,, i.e. with viscosity much larger 
than resistivity, are widely spread. In particular, the proper- 
ties of turbulence in high Pr„, fluid are relevant both to fully 
and partially ionized media of the Local Bubble. 

For high Pr,„ fluid Cho, Lazarian & Vishniac (2002b) re- 
ported a new regime of MHD turbulence. Lazarian, Vishniac 
& Cho (2004) showed that while the spectrum of volume- 
averaged magnetic fluctuations scales as E[,{k) ^ k^^ (see 
FigEJ, the pressure within intermittent magnetic structures 
increases with the decrease of the scale k and the filling 
factor decreases oc k^^. The magnetic pressure compresses 
the gas as demonstrated in Fig. |3] More importantly, ex- 
tended current sheets that naturally emerge as magnetic field 

^ We, however, warn our reader not to identify the viscosity and the 
effects of neutral gas on turbulence. The actual physics of neutral-ion 
interactions in turbulent gas goes beyond the effects of inducing the 
drag. 
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fluctuates in the plane perpendicular to the mean magnetic 
field (see Fig.O. It was speculated in Lazarian (2007a) that 
these current sheets can account for the origin of the small 
ionized and neutral structures (SINS) on AU spatial scales 
(Dieter et al. 1976, Heiles 1997, Stanimirovic et al. 2004). 

Goldreich & Sridhar (2006) appealed not to high Pr„, 
MHD turbulence per se, but to the generation of the mag- 
netic field in the turbulent plasma (see Schekochihin et al. 
2004) to account for the high amplitude, but small scale 
fluctuations of plasma density observed in the direction of 
the Galactic center. We believe that the regime of dynamo in 
Schekochihin et al (2004) and the turbulence in Lazarian et 
al. (2004) have similarities in terms of the density enhance- 
ment that are created. Although in the case of magnetic tur- 
bulence with sufficiently strong mean magnetic field, global 
reversals, that Goldreich & Shrindhar (2006) appeal to in 
compressing plasma, do not happen, the reversals of the mag- 
netic field direction occur in the direction perpendicular to 
the mean magnetic field. As the mean magnetic field goes 
to zero, the two regimes get indistinguishable. A systematic 
study of the density enhancement within the Local Bubble, 
where more detailed knowledge of the plasma and magnetic 
field properties are potentially available, may test the theo- 
retical constructions above. 



3 Intermittency: small volumes with extreme conditions 

An anisotropic spectrum alone, say £'(k)iik, cannot charac- 
terize MHD turbulence in all its complexity because it in- 
volves only the averaged energy in motions along a particu- 
lar direction. To have a full statistical description, one needs 
to know not only the averaged spectrum of a physical vari- 
able but higher orders as well. The tendency of fluctuations 
to become relatively more violent but increasingly sparse in 
time and space as the scales get smaller, so that their influ- 
ence remains appreciable, is called intermittency . The inter- 
mittency increases with the ratio of the size scales of injec- 
tion and dissipation of energy, so the very limited range of 
scales within numerical simulations may fail to reflect the 
actual small scale processes. 

Falgarone et al. (2005, 2006, 2007) and collaborators 
(Hily-Blant & Falgarone 2007, Hily-Blant et al. 2007 and 
refs. therein) attracted the attention of the interstellar com- 
munity to the potential important implications of intermit- 
tency. A small and transient volume with high temperatures 
or violent turbulence can have significant effects on the net 
rates of processes within the ISM. For instance, many in- 
terstellar chemical reactions (e.g., the strongly endothermic 
formation of CH+) might take place within very intensive 
intermittent vortices. The aforementioned authors claimed 
the existence of the observational evidence for such reac- 
tions and heating. If the effects of intemittency are as strong 



as Falgarone et al. believe, they should not be neglected 
when processes in the Local Bubble are considered. 

To get insight into intemittency effects, it is necessary 
to study high moments of velocity fluctuations. Fortunately, 
both laboratory and numerical studies demonstrate that the 
higher moments of velocity fluctuations can be predicted re- 
markably well by the expression derived by She and Lev- 
eque (1994) (see discussion in Lazarian 2006a and ref. therein) 
The key parameter is the spatial dimension of the dissipa- 
tion structures, D, ^ \ for filamentary vortices and ^ 2 for 
sheets in MHD (see MuUer & Biskamp 2000, Cho, Lazarian 
& Vishniac 2002, Boldyrev et al. 2002). In some instances 
D can be a fractal dimension in between. The dimension can 
be obtained even with low resolution numerics. This opens 
an avenue of evaluating the effects of intermittencies, such 
as intense local heating or transient effects, for Re and Rm 
numbers that will probably never be achieved in numerical 
simulations. 

The She-Leveque (1994) model describes the p\h power 
of longitudinal velocity fluctuations at scale I, i.e. (SV/)'' ~ 
l^p, where 5Vi = 5V/1/Z. Here 5Vi is the mean difference 
of the local velocity from the average. For non-intermittent 
Kolmogorov turbulence, — p/3 (recall the well-known 
result, 5Vi ^ 1^1^). In the She-Leveque (1994) model of tur- 
bulence, is a more complex function that depends on 
(a) the scaHng of velocity 5Vi ^ I", (b) the energy cascad- 
ing rate f/ ^ 1^1^, and (c) the dimension of the dissipation 
structures, D. The She-Leveque (1994) expression is ^p = 
ai3(l-j3) + (3-D)(l-[l-j3/(3-D)]«'^),witha-l/3 
and /3 ^ 2/3 for MHD turbulence. The principal dependence 
of £,p is on D, between 1 and 2. 

DubruUe (1994), She & Waymire (1995) have figured 
out that the scaling exponents ^p of the She-Leveque model 
correspond to so-called generalized log-Poisson distribution 
of the local dissipation rate £/. With this input, if we con- 
strain the model parameters from low resolution numerical 
simulations we can calculate the probability of a fluctuation 
exceeding a given threshold of the local rate of deposition of 
energy, for arbitrary high values of Re and Rm. Very impor- 
tantly, predictions can be made for the real ISM! 

Fig. |4] shows how to determine D. It involves a simu- 
lation (Mach numbers shown) for the component of V per- 
pendicular and also parallel to the local B. The dot-dashed 
and dashed lines are the respective predictions of the She- 
Leveque (1994) expression, with the values of D as indi- 
cated. The simulation clearly provides appropriate values of 
D for the two components. 

The bottom part of Figure|4]shows calculations in Beres- 
nyak & Lazarian (2008c) for the volume fractions of vari- 
ous dissipation rates (i.e., heating). While the temperatures 
achieved will depend upon the cooling functions, some im- 
portant conclusions are available from the analysis of Figure 
m Indeed, the model of chemical reactions by Falgarone et 
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Fig. 3 Left: Filaments of density created by magnetic compression of tlie gas in tlie viscosity-damped regime of MHD turbulence. Darker regions 
correspond to higher density. The viscous damping scale l,. is much larger than the cuiTcnt sheet thickness d. This creates large observed density 
contrasts. Center. Magnetic reversals (in the plane ± to mean (B)) that create compressions of density. Darker regions correspond to higher 
magnetic field. Right: Spectra of density and magnetic field are similar, while velocity is damped. The resistive scale in this regime is not L/Rm 
but LRnC^I^. From Beresnyak & Lazarian. in prep.) 
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Fig. 4 Left Panel: The intermittencies of velocities in our subAlfvenic, 
Ma = 0.7 supersonic Ms = 7 MHD simulations. From Kowal & Lazar- 
ian (2008). Right Panel: Volume fraction with the dissipation rate is 
higher than the mean rate for the She-Leveque model of intermittency 
with D = 1 and 2. From Beresnyak & Lazarian. 

al. requires that a substantial part of the turbulent cascade 
energy dissipate in the very intermittent structures. Figure|4] 
shows that the bulk of the energy dissipates within structures 
where the dissipation rate is higher than the mean value less 
than the factor of 100, provided that the She-Leveque model 
is valid. This provides stringent constraints on what chem- 
istry we could expect to be induced by intermittent turbulent 
heating. 

Interestingly enough, the case of intermittency studies 
supports our point of the futility of the "brute force" numer- 
ical approach. For instance, for a typical ISM injection scale 



of 50 pc, the Reynolds number can be as high as Re = 10 . 
In comparison, numerical simulations provide ^ Re^ boxes 
for the present record resolution of a hydrodynamic simula- 
tion with 4096^* boxes. 

Even if the intermittency of turbulence is not as impor- 
tant as is in Falgarone et al. model, its implications for both 
ISM and Local Bubble plasma may be very important. We 
see in Figure |4] that 10% of the energy deposition is local- 
ized in just 10^^% of the volume! Such a concentration of 
energy dissipation in high Re can have many important con- 
sequences on inhomogeneous heating (see Sonnentrucker et 
al. 2006). It is clear that further studies in this direction are 
necessary. 



4 Modification of turbulence by cosmic rays 

It is easy to argue that the interaction of turbulence and cos- 
mic rays (CRs) is of the most important processes where 
properties of MHD turbulence are essential. Indeed, the in- 
teraction of turbulence with CRs is a cornerstone of CRs 
propagation and acceleration models (e.g., Ginzburg 1966, 
Jokipii 1966, Wentzel 1969, Schlikeiser 2002 and ref. therein). 
To account for the interaction properly, one must know both 
the scaling of turbulence and the interactions of turbulence 
with various waves produced by CRs. As we mentioned above, 
slab Alfvenic modes, which are an essential part of many 
models of cosmic ray propagation (see Goldstein 1976, Bieber 
et al. 1994, Matthaeus et al. 2003, Shalchi, Bieber & Matthaeus 
2004, Shalchi 2005, 2006 and references therein), are not 
observed in direct MHD simulations. This does not mean 
that these modes are necessarily absent in realistic astro- 
physical environments. For instance, the influence of the 
CRs back onto the turbulence may be important. 

The problem is that energetic particles or CRs being an 
important component in both the ISM and Local Bubble are 
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not a part of the simulations of MHD turbulence. Neverthe- 
less, CRs are dynamically important and well coupled to the 
rest of the ISM through magnetic fields and magnetic tur- 
bulence. The spectra of both CRs and interstellar turbulence 
show nice power laws, leading R. Jokipii (2001) to suggest a 
strong interrelation of the two. It has to be noted that the to- 
tal pressure of CRs in the ISM is of the order of the kinetic 
pressure pv^ and, in many cases, exceed the thermal pres- 
sure by as much as a factor of 10. One cannot exclude that 
CRs may be as important for ISM simulations as dark matter 
is for the problems of galaxy formation and dynamics. 

The direct application of the results of MHD turbulence 
theory to the problems of the CRs propagation and accel- 
eration is well justified only when the modification of tur- 
bulence by CRs can be neglected. This has caused contro- 
versies since the classical work by Parker (1965). CRs can 
induce and are influenced by two instabilities: (a) streaming 
instability (Wentzel 1969, 1974), widely discussed in the CR 
literature and frequently referred to as the source of the pos- 
tulated scattering "slab" , i.e. with k||B Alfvenic component, 
and (b) the anisotropic pressure kinetic gyroresonance in- 
stabiHty (see Mikhailovskii 1975, Gary 1993, Kulsrud 2004 
and refs. therein), well known to the plasma community. 

The streaming instability, as implied by its name, re- 
quires an anisotropic distribution of particles, while in the 
interstellar gas, if we judge from our local measured values, 
the distribution of CRs is close to isotropic. In addition, Yan 
& Lazarian (2002) noticed that the streaming instability can 
be non-linearly suppressed by ambient turbulence. Farmer & 
Goldreich (2004) estimated that galactic turbulence can effi- 
ciently suppress the formation of the "slab" Alfvenic waves 
for protons with relativistic factor 7 > 100. Numerical cal- 
culations in Beresnyak & Lazarian (2008a) confirmed (see 
Fig EJ those predictions. 

Lazarian & Beresnyak (2006) research showed that com- 
pressible MHD turbulence can generate kinetic gyroreso- 
nance instability that can couple efficiently CRs with tur- 
bulence and backreact on turbulence, modifying its spec- 
trum (see Fig. |6]l. The instability arises from compressing 
magnetic fields with gyrating CRs. Indeed, the CR pressure 
becomes anisotropic as a consequence of adiabatic invari- 
ant conservation, thereby inducing the instability (Lazarian 
& Beresnyak 2006). The latter gives rise to "slab" Alfvenic 
perturbations that scatter and randomize the CR momenta 
via gyroresonance. As the "slab" mode grows, it increases 
CR scattering and coupling of CRs and magnetic fields. 

Further research of the kinetic gyroresonance instabili- 
ties, as well as possibly other instabilities induced by CRs, 
should improve our understanding of the importance of the 
CRs-generated slab Alfven waves. Unlike earlier models which 
used somewhat arbitrary (see §2) prescriptions for the am- 
plitude of the "slab" Alfvenic mode (see Bieber et al. 1994), 
future theories appealing to such modes should be able to 
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Fig. 5 Decorrelation of a plane, = Alfven wave by turbulence. 
Lower picture shows the energy density of a wave in cylindrical k- 
space. Alfven waves were injected at ^|| = 17. From Beresnyak & 
Lazarian (2008a). 
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Fig. 6 Energy density of compressive modes and Alfvenic slab-type 
waves, induced by CRs. The energy is transferred from the mean free 
path scale to the CR Larmor radius scale. See more on feedback in 
Lazarian & Beresnyak (2006). 

evaluate their expected amplitudes (see first attempts to do 
this in Lazarian & Beresnyak 2006). We should also men- 
tion, that the slab mode will be subject to damping while 
interacting with the surrounding strong Alfvenic turbulence 
(Yan & Lazarian 2002, Farmer & Goldreich 2004, Lazarian 
& Beresnyak 2006, Beresnyak & Lazarian 2008). Includ- 
ing this damping in Shalchi, Lazarian & Schlickeiser (2007) 
improved the fit to the solar wind measurement by Ulysses 
(Gloeckler et al. 1995) and AMPTE spacecraft (Mobius et 
al. 1998). 

5 Selected implications of magnetic turbulence 

The implications of magnetic turbulence for astrophysical 
fluids, e.g. interstellar medium, are numerous. For instance, 
turbulence can heat Diffuse Ionized Gas within the Milky 
Way (see Minter & Spangler 1997, Cho et al. 2002), deter- 
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mine the evolution of molecular clouds (see McKee & Os- 
triker 2007), scatter and accelerate cosmic rays (see Schlick- 
eiser 2002). Below we consider a couple of selected exam- 
ples relevant to the research by the authors to illustrate the 
progress and problems on the way of evaluating the effects 
of turbulence. 



5.1 Cosmic ray scattering by fast modes of MHD 
turbulence 

As we mentioned earlier, numerical simulations of MHD 
turbulence supported the GS95 model of turbulence, which 
does not have the "slab" Alfvenic modes that produced most 
of the scattering in the earlier models of CR propagation. 
Can the turbulence that does not appeal to CRs back-reaction 
(see §4) produce efficient scattering? 

In the models of ISM turbulence (Armstrong et al. 1994, 
McKee & Ostriker 2007), where the injection happens at 
large scale, fast modes were identified as a scattering agent 
for cosmic rays in interstellar medium Yan & Lazarian (2002, 
2004). These works made use of the quantitative description 
of turbulence obtained in CL02 to calculate the scattering 
rate of cosmic rays. The results are shown in Fig.|2] For in- 
stance, the scattering rate of relativistic protons by Alfvenic 
turbulence was shown to be nearly lO'" times lower than 
the generally accepted estimates obtained assuming the Kol- 
mogorov scaling of turbulence. Although this estimate is lO'* 
times larger than that obtained by Chandran (2000), who 
employed GS95 ideas of anisotropy, but lacked the quanti- 
tative description of the eddies, it is clear that for most inter- 
stellar circumstances the Alfvenic scattering is suppressed. 
The low efficiency of scattering by Alfvenic modes arise 
from high anisotropy of the modes at the scales of cosmic 
ray gyroradius. 

YL02 showed that the scattering by fast modes, which 
are isotropic (CL02), dominates (see Fig. |7]i. However, fast 
modes are subject to both collisional and coUisionless damp- 
ing which was taken into account in Yan & Lazarian (2004). 
More recent studies of cosmic ray propagation and acceler- 
ation that explicitly appeal to the effect of the fast modes in- 
clude Cassano & Brunetti 2005, Brunetti 2006, Brunetti & 
Lazarian (2007), Yan & Lazarian (2008) and Yan, Lazarian 
& Petrosian (2008). Incidentally, fast modes have been also 
identified as primary agents for the acceleration of charged 
dust particles (Yan & Lazarian 2003, Yan, Lazarian & Draine 
2004). 

^ On the basis of weak turbulence theory, Chandran (2005) has ar- 
gued that high-frequency fast waves, which move mostly parallel to 
magnetic field, generate Alfven waves also moving mostly parallel to 
magnetic field. We expect that the scattering by thus generated Alfven 
modes to be similar to the scattering by the fast modes created by them. 
Therefore we expect that the simplified approach adopted in Yan & 
Lazarian (2004) and the papers that followed to hold. 



5.2 Turbulent reconnection and cosmic ray acceleration 

Magnetic reconnection can be associated with the ability of 
magnetic flux tubes to change their topology, while being 
submerged within conducting fluids (see Biskamp 1996). 
Due to high numerical diffusivity of present-day simula- 
tions, reconnection is fast there, which, for instance, means 
that magnetic fields in the ISM, Local Bubble and solar wind 
simulations change their topology fast. Is this true for the 
real astrophysical circumstances? 

Recent years have been marked by a substantial progress 
in simulations of coUisionless reconnection (see Shay & Drake 
1998, Bhattacharjee et al. 2003, Drake et al. 2006 and refer- 
ences therein). However, while the researchers argue whether 
Hall MHD or fully kinetic description (Daughton 2006) is 
necessary, one statement is definitely true. If magnetic re- 
connection is only fast in coUisionless environments, most 
of the MHD simulations, e.g. of interstellar medium, accre- 
tion disks, stars, where the environment is collisional, are in 
error. We shall argue below that this radical conclusion may 
not be true and the reconnection is also fast in most astro- 
physical collisional environments. 

Lazarian & Vishniac (1999, henceforth LV99) consid- 
ered turbulence as the agent that makes magnetic recon- 
nection fast. The scheme proposed in LV99 there differs 
appreciably from the earlier attempts to enhance reconnec- 
tion via turbulence (Speiser 1970, Jacobson & Moses 1984, 
Matthaeus & Lamkin 1985, Bhattacharjee & Hameiri 1986, 
Hameiri & Bhattacharjee 1987, Straus 1988, see Lazarian et 
al. 2004 for a detailed comparison). The scheme proposed 
is a generalization of the Sweet-Parker scheme (see Fig. [8]). 
The problem of the Sweet-Parker model is that the reconnec- 
tion is negligibly slow for any realistic astrophysical con- 
ditions. However, astrophysical magnetic fields are generi- 
cally turbulent. 

LV99 consider the case in which there exists a large 
scale, well-ordered magnetic field, of the kind that is nor- 
mally used as a starting point for discussions of reconnec- 
tion. In addition, we expect that the field has some small 
scale 'wandering' of the field lines. On any given scale the 
typical angle by which field lines differ from their neighbors 
is <C 1, and this angle persists for a distance along the field 
lines A|| with a correlation distance X± across field lines (see 
Fig.©. 

The modification of the global constraint induced by mass 
conservation in the presence of a stochastic magnetic field 
component is self-evident. Instead of being squeezed from a 
layer whose width is determined by Ohmic diffusion, the 
plasma may diffuse through a much broader layer, Ly ~ 
(^2^/2 (ggg Pig [8]|^ determined by the diffusion of magnetic 
field lines. This suggests an upper limit on the reconnection 
speed of ^ VA((y^)'/^/Lv). This will be the actual speed of 
reconnection if the progress of reconnection in the current 
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Fig. 7 Implications of Interstellar Turbulence for the Cosmic Rays and Interstellar Dust. Upper left: Rate of CR scattering by Alfven waves 
versus CR energy. The lines at the top of the figure are the accepted estimates obtained for Kolmogorov turbulence. The dotted curve is from 
Chandran (2000). The analytical calculations are given by the solid line with our numerical calculations given by crosses. Upper right: The rate 
of CR scattering (v) by fast modes in magnetically dominated plasma. The rate of scattering depends on damping of the fast waves (see , which 
in turn depends on the ratio of gaseous to magnetic pressure (j8 = Pgus/Pmag)- Lower left: Individual trajectories of CRs tracked by the Monte 
Carlo scattering code. B is obtained through 3-D simulations of MHD turbulence. These calculations provide estimates of CR diffusion. Lower 
right: Velocities of charged dust grains in cold neutral media (CNM). Gyroresonance with fast modes ("fast gyro") dominates for large grains, 
while hydrodynamic drag ("fast hydro") for small grains. The cutoffs for Alfven and fast gyro (vertical lines) are due to MHD turbulence damping 
caused by neutral-ion collisions. From Yan & Lazarian 2002, 2003, 2008. 



sheet does not impose a smaller limit. The value of (y^)'/^ 
can be determined once a particular model of turbulence is 
adopted, but it is obvious from the very beginning that this 
value is determined by field wandering rather than Ohmic 
diffusion as in the Sweet-Parker case. 

What about limits on the speed of reconnection that arise 
from considering the structure of the current sheet? In the 
presence of a stochastic field component, magnetic recon- 
nection dissipates field lines not over their entire length ~ Lx 
but only over a scale A|| <C Lx (see Fig. [SJ, which is the 
scale over which magnetic field line deviates from its orig- 
inal direction by the thickness of the Ohmic diffusion layer 
Xj^ « r\ /Vrec.iocai- If the angle of field deviation does not 
depend on the scale, the local reconnection velocity would 
be ^Va^ and would not depend on resistivity. In LV99 it 
is taken into account that does depend on scale (see §2). 
Therefore the local reconnection rate Vrec.iocai is given by 
the usual Sweet-Parker formula but with A|| instead of Lx, 
i.e. Vrec.iocai ~ Va{VaX\\/ 7])^^l^ . It is obvious from Fig. [8] 
that ^ Lx/X\\ magnetic field lines will undergo reconnection 
simultaneously (compared to a one by one line reconnec- 



tion process for the Sweet-Parker scheme). Therefore the 
overall reconnection rate may be as large as Vrec,giobai ~ 
V'A(i'.r/A||)(V'AA||/?7)^'/^. Whether or not this limit is impor- 
tant depends on the value of A|| . 

The relevant values of A|| and (y^)^^^ depend on the mag- 
netic field statistics. This calculation was performed in LV99 
using the GS95 model of MHD turbulence providing the up- 
per limit on the reconnection speed: 



Vr.up Va min 



(1) 



where I and v/ are the energy injection scale and turbulent 
velocity at this scale respectively. In LV99 other processes 
that can impede reconnection were found to be less restric- 
tive. For instance, the tangle of reconnection field lines cross- 
ing the current sheet will need to reconnect repeatedly be- 
fore individual flux elements can leave the current sheet be- 
hind. The rate at which this occurs can be estimated by as- 
suming that it constitutes the real bottleneck in reconnection 
events, and then analyzing each flux element reconnection 
as part of a self-similar system of such events. This turns out 
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Fig. 8 Upper plot: Sweet-Parker model of reconnection. The outflow 
is limited by a thin slot A , which is determined by Ohmic dilTusivity. 
The other scale is an astrophysical scale A. Middle plot: Turbu- 
lent reconnection model that accounts for the stochasticity of magnetic 
field lines. The outflow is limited by the diffusion of magnetic field 
lines, which depends on field line stochasticity. Low plot: An individual 
small scale reconnection region. The reconnection over small patches 
of magnetic field determines the local reconnection rate. The global 
reconnection rate is substantially larger as many independent patches 
come together. 



to impede the reconnection. As the result, LV99 concludes 
that ^ is not only an upper limit, but is the best estimate 
of the speed of reconnection. The model has been recently 
tested numerically in Kowal et al. (2008) (see Fig. |9]l. The 
thick current sheets observed during the 2003 November 4 
Coronal Mass Ejection reported in Ciaravella & Raymond 
(2008) are also consistent with the LV99 model. 

The most interesting process is the first-order Fermi ac- 
celeration that is intrinsic to the turbulent reconnection. To 
understand it consider a particle entrained on a reconnected 
magnetic field line (see Figur^Sll. This particle may bounce 
back and forth between magnetic mirrors formed by oppo- 
sitely directed magnetic fluxes moving towards each other 
with the velocity Vr. Each of such bouncing will increase 
the energy of a particle in a way consistent with the require- 
ments of the first-order Fermi process. The interesting prop- 
erty of this mechanism that potentially can be used to test 
the idea observationally is that the resulting spectrum is dif- 
ferent from those arising from shocks. Gouveia Dal Pino & 
Lazarian (2003) used this mechanism of particle accelera- 



tiorOto explain the synchrotron power-law spectrum arising 
from the flares of the microquasar GRS 1915H-105. Note, 
that the mechanism acts in the Sweet-Parker scheme as well 
as in the scheme of turbulent reconnection. However, in the 
former the rates of reconnection and therefore the efficiency 
of acceleration are marginal in most cases. 

5.3 Perpendicular diffusion and subdiffusion 

Modeling of clouds submerged in hot plasma present one of 
the challenges to the research of the Local Bubble physics. 
It is easy to see that the diffusion heat perpendicular to mag- 
netic field is very slow if the field is laminar Incidentally, 
the same is true in terms of CR diffusion, while the obser- 
vations indicate that the diffusion of the CRs perpendicular 
to magnetic field in the Galaxy is reduced just by a factor 
~ 3 compared to the parallel diffusion (Jokipii 1999). How 
could this be? 

The reader may already guess that, similar to the re- 
connection problem discussed above, magnetic field wan- 
dering may allow the particles to diffuse (Jokipii & Parker 
1969, Giacalone, J., & Jokipii, J. R. 1999). Analogously 
with the reconnection example, one may expect that the ex- 
pression for the field wandering would depend on the model 
of turbulence accepted. Within the GS95 model of turbu- 
lence and assuming that the injection turbulent velocity is 
equal to the Alfven one, the calculations for heat diffusion 
were performed in Narayan & Medvedev 2001 A general 
case of arbitrary injection velocity was considered in Lazar- 
ian (2006). The heat diffusion in plasma was found to be 
a function of the Alfven Mach number defined as the 
ratio of the turbulent velocity at the injection scale to the 
Alfven speed in the medium. For the case of strong mean 
magnetic field, i.e. forM^ < 1, and mean free path of a elec- 
tron X less than the injection scale of the turbulence L, the 
heat diffusion coefficient was obtained to be 1 /'iM'^KspUzer^ 
where K^pitzer — ■^I'r/ie™ is the usual Spitzer heat diffusion 
coefficient of unmagnetized plasma. The factor 1 /3 there re- 
flects the ID nature of diffusion along magnetic field lines, 
while the poweJ*^ reflects the inefficiency of magnetic 
field line wandering in the presence of the strong mean mag- 
netic field. Similarly, for superAlfvenic turbulence, i.e. for 
Ma > 1 one gets the heat diffusion coefficient l/3Kspitzer 
if the mean free path of an electron A < LM^^ and gets 
1 /3{LM^^ / X)Kspitzer in the opposite regime. 

^ The mechanism has physical similarities to the acceleration mech- 
anism that was proposed later for electrons by Drake et al. (2006b). In 
Drake's mechanism, similarly, to the Matthaeus, Ambrosiano & Gold- 
stein (1984) mechanism, however, the process of acceleration hap- 
pens within 2D contracting loops. For LV99 model of reconnection 
the generic configuration of magnetic field are contracting spirals. 

There is a wrong power in the original Lazarian (2006) paper, 
which was corrected in Lazarian (2007). 



13 




Fig. 9 Left panel Testing of the model of 3D turbulent reconnection with MHD simulations in Kowal et al. (2008) The results, consistently with 
the LV99 model, show an increase of the reconnection rate with the increase of the injection scale of the turbulence. The reconnection is bursty, as 
is expected for the turbulence model. The reconnection rate is normalized to the Alfvenic velocity, time is given in the units of Alfvenic crossing 
times. Initially the turbulence is slow, but its velocity increases as the turbulence is injected. Central Panel. The same as for the left panel, but 
the injection energy is being changed. The results are also consistent with the LV99 model. Right panel: Cosmic rays spiral about a reconnected 
magnetic field line and bounce back at points A and B . The reconnected regions move towards each other with the reconnection velocity Vr . The 
advection of cosmic rays entrained on magnetic field lines happens at the outflow velocity, which is in most cases of the order of Va. Bouncing at 
points A and B happens because either of streaming instability or turbulence in the reconnection region. 



Cho et al. (2003) noticed that in magnetized turbulence, 
the turbulence eddies should induce in plasma the heat ad- 
vection with the effective diffusion coefficient 2 /SLVl, where 
the factor of 1 /3 is similar to the advective diffusion in hy- 
drodynamics (Lesieur 1990) and the factor of 2 takes into ac- 
count that both electrons and protons participate in the pro- 
cess. The relative importance of the two diffusion processes 
is exemplified in Figure [Tol where the parameter space for 
the dominance of the heat advection by turbulent motions 
and the dominance of the heat transfer by electrons are de- 
fined. These results for the fully ionized plasma can be used 
for parameterizing heat transport in the astrophysical codes. 
They also can easily be generalized for the partially ionized 
gas. Note, that following individual particles in simulations 
may be prohibitively expensive (see Tilley & Balsara 2006). 
Moreover, the actual structure of magnetic field is distorted 
by the limited numerical resolution. 

The last point of magnetic field structure in numerical 
simulations is also important for the CRs diffusion perpen- 
dicular to magnetic field. Take, for instance, the issue of sub- 
diffusion, which is the non-diffusive behavior arising from 
the CRs retracing their trajectories, as a result of the back- 
ward scattering (see Kota & Jokipii 2000, Qin et al. 2002, 
Webb et al. 2006). In the case of the unrealistic turbulence 
with only one scale of turbulent motions Iturhi the retrac- 
ing of particles stops as soon as the particles diffuse the 
so-called Rechester & Rosenbluth (1978) length, which is 
kurb ^^{kurb I fUir), where ri^,,- is the Larmor radius of a charged 
particle (see also Chandran & Cowley 1998). In the case of 
realistic turbulence with a range of scales one should use the 
dissipation scale instead of the f,„ri (Narayan & Medvedev 
2001, Lazarian 2006, Yan & Lazarian 2008). The correspond- 
ing scale for the Alfvenic turbulence in the ISM may be less 
than Icrit ~ 10^ km, thus one should expect subdiffusion only 
at scales less than this. For the solar wind turbulence Icrk is 
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Fig. 10 The parameter spaces for the advection of heat by turbulent 
motions and heat conduction by plasma electrons. Sonic Mach number 
M, is plotted against the Alfven Mach number . The mean free path 
of electrons is 1, while a is the square root of the electron to proton 
mass ratio, i.e. (m^/mpYI^, j3 is a numerical factor fa 4. From Lazarian 
2006, 2007. 



even smaller, of the order of 10^ km. On the contrary, for 
the numerical simulations with a limited inertial range, the 
subdiffusion may be much more important (see Qin et al. 
2002). 



6 Studying magnetic turbulence in Interplanetary 
medium using comets 

As we illustrated with selected examples above, the model- 
ing of turbulent astrophysical fluids is far from being simple. 
Therefore, observational testing is essential. We are fortu- 
nate to have in situ measurements of magnetic field in the 
interplanetary medium (see Opher et al. 2007). The advan- 
tage of direct studies of magnetic perturbations by space- 
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crafts has been explored through many important missions. 
Such studies, unHke numerical ones, may deliver informa- 
tion about the actual magnetic turbulence at high Re and Rm 
numbers. However, the spacecraft measurement are rather 
expensive. Are there any other cost-effective ways to study 
magnetic turbulence in interplanetary medium? 

Below we present a new way to explore the turbulence 
in interplanetary medium by using the alignment of sodium 
atoms ejected from comets. Atomic alignment of atoms in 
their ground state was an effect studied in the middle of 
the previous century in laboratories, in relation to maser 
research (Hawkins 1955, Kastler 1957). The alignment of 
atoms is understood in terms of their non-equilibrium distri- 
bution of ground level substates. Therefore, to be aligned 
in the ground state the atom in question should have ei- 
ther fine or hyper-fine structure. The non-equilibrium dis- 
tribution arises from radiative pumping, while the magnetic 
field realigns atoms inducing their Larmor precession. As 
the substates of the ground level are long-lived, even a small 
magnetic field modifies the distribution and therefore the po- 
larization arising from the aligned atoms (see Yan & Lazar- 
ian 2007 for more details). 

Comets are known to have Sodium tails and Sodium is 
an atom that can be aligned by radiation and realigned by 
solar wind magnetic fields. This opens an opportunity of 
studying magnetic fields in the solar wind from the ground, 
but tracing the polarization of the Sodium line. At the mo- 
ment this is a suggestion supported by the synthetic ground 
based observations. For our synthetic observations we used a 
model of magnetic field in Liu et al. (2008). It employed the 
space weather model developed by University of Michigan, 
namely. The Space Weather Modeling Framework (SWMG) 
(Toth et al 2005). More specifically, this is a solar corona 
and inner heliosphere model that extends the description of 
media from the solar surface to 1 AU. 

The structure of the magnetic field in the heliosphere can 
be studied by the polarization of Sodium D2 emission in the 
comet's wake. Though the abundance of sodium in comets 
is very low, its high efficiency in scattering sunlight makes it 
a good tracer (Thomas 1992). As discussed in Yan & Lazar- 
ian (2007), the gaseous sodium atoms in the comet's tail ac- 
quire angular momentum from the solar radiation, i.e. they 
are aligned. Resonant scattering from these aligned atoms 
is polarized. Distant from comets, the Sun can be consid- 
ered a point source. As shown in Fig. (TT] the geometry of 
the scattering is well defined, i.e., the scattering angle Oq is 
known. The alignment is modulated by the local magnetic 
field. The polarization of the sodium emission thus provides 
exclusive information on the magnetic field in the interplan- 
etary medium. We take the data cube from the spacecraft 
measurement as described above. Depending on its direc- 
tion, the embedded magnetic field alters the degree of align- 
ment and therefore the polarization of the light scattered by 



the aligned atoms. Fig[TT|illustrates the trajectory of a comet 
along which the magnetic field varies and the polarization 
of Sodium D2 emission changes accordingly. By comparing 
observations with it, we can determine the structure of mag- 
netic field in the heliosphere. For interplanetary studies, one 
can investigate not only spatial, but also temporal variations 
of magnetic fields. Since alignment happens at a time scale 
Tr, magnetic field variations on this time scale will be re- 
flected. This can allow for a cost-effective way of studying 
interplanetary magnetic turbulence at different scales. On 
the basis of the results above we expect that comets may 
become an important source of information about interplan- 
etary magnetic fields and their variations. 



7 Discussion and Summary 

The goals of the review above were, first of all, to demon- 
strate the necessity of detailed studies of magnetic turbu- 
lence for modeling of various astrophysical environments, 
in particular, the environment of the Local Bubble and, sec- 
ondly, to appeal for the necessity of new ways of studying 
magnetic fields in solar wind. To address the former goal 
we discussed the properties of MHD turbulence in different 
environments (fully ionized and partially ionized gas) for 
different ways of excitation (balanced and imbalanced tur- 
bulence). The second goal was addressed by considering a 
new promising way of tracing magnetic field structure in the 
interplanetary medium. 

The examples are representative of authors' interests, but 
they vividly show that the "brute force approach" to sim- 
ulating Local Bubble and other astrophysical environments 
may easily fail. For instance, we showed that it would not be 
fruitful for calculating of the intermittency of environmen- 
tal heating and it is likely to fail in representing the effects 
of subdiffusion. It also may deliver unreliable results when 
simulating heat transfer in magnetized plasmas etc. There- 
fore detailed quantitative modeling may require, first of all, 
the creation of the "tool box" of the particular recepiees of 
how to parameterize particular properties of turbulent fluid. 
We claimed that studies of these properties were not only de- 
sirable, but were absolutely necessary. For instance, unless 
we reach a consensus on the rate of reconnection in colli- 
sional environments, the simulations of the interstellar gas 
will be highly suspect (see the corresponding discussion in 
§5.2). All in all, we believe that the progress can be achieved 
via better understanding of fundamental properties of mag- 
netized plasma and aggressive testing the results of mod- 
eling with observations. The fact that the properties of tur- 
bulent fluid are important for a wide range of astrophysical 
problems makes it easier to concentrate efforts and resources 
on such studies. This also calls for more vigorous scientific 
exchanges between the disciplines. 
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Fig. 11 Left Panel. Schematic of the resonance scattering of sunlight by the sodium in comet wake. The sodium tail is points in the direction 
opposite to the Sun. The observer on the Earth sees the stream at the angle Oq. Magnetic field realigns atoms via fast Larmor precession. Thus the 
polarization traces the interplanetary magnetic fields. Right Panel. The magnetic field from the model of inner heliosphere (blue lines) (Toth et al. 
2005) and the predicted vectors of linear polarization of Sodium D2 emission (red lines) along the comet trajactory. From Yan et al. 2008. 



Our points of the review above can be briefly summa- 
rized in the foflowing way: 

1. Properties of astrophysical fluids, including those of so- 
lar wind, Local Bubble and interstellar turbulence are, to 
large extend, determined by magnetic turbulence. The "brute 
force" approach aimed at detailed modeling of the afore- 
mentioned astrophysical environments keeping magnetic tur- 
bulence realistic is doomed. Instead, we feel fliat focused 
studies of particular physical phenomena in turbulent mag- 
netized plasma can clarify when the simulations that take 
into account many processes at once provide a correct phys- 
ical picture. These studies can help parameterize effects of 
turbulence in numerical codes. 

2. By considering imbalanced turbulence, as this is the case 
of the Solar wind turbulence, and high-Pr number turbu- 
lence, that can approximate the turbulence in the partially 
ionized plasma, we demonstrated that magnetic turbulence 
can have many different properties. These properties affect 
the transport properties of the magnetized media, e.g. scat- 
tering of cosmic rays, the formation of density enhance- 
ments, including the formation of small ionized and neu- 
tral structures (SINS), the change of magnetic topology, i.e. 
magnetic reconnection, thermal conductivity. 

3. We advocate the approach to simulating processes in the 
aforementioned environments that uses the synergy of fo- 
cused studies of particular processes in turbulent magnetized 
plasmas, numerical simulations of astrophysical situations 
and observational studies. We also feel that collecting of 
new data on magnetic fields is essential. In the current situ- 
ation, it is promising to explore new techniques of magnetic 
field studies, e.g. the ones that makes use of atoms aligned in 
their ground state by radiation and realigned by the external 
magnetic field. In particular, we discuss possible studies of 
interplanetary turbulence using Na aligned atoms in a comet 
wake. 
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